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1. Summary 

As part of the work performed under NASA contract # NAS5-32648, we have computed the 

3- point and 4-point correlation functions of the COBE-DMR 2-year and 4-yea.r anisotropy maps. 
The motivation for this study was to search for evidence of non-Gaussian statistical fluctuations 
in the temperature maps: skewness or asymmetry in the case of the 3-point function, kurtosis 
in the case of the 4-point function. Such behavior would have very significant implications for 
our understanding of the processes of galaxy formation, because our current models of galaxy 
formation predict that non-Gaussian features should not be present in the DMR maps. The results 
of our work showed that the 3-point correlation function is consistent with zero and that the 

4- point function is not a very sensitive probe of non-Gaussian behavior in the COBE-DMR data. 

Our computation and analysis of 3-point correlations in the 2-year DMR maps was published 
in the Astrophysical Journal Letters, volume 446, page L67, 1995. Our computation and analysis 
of 3-point correlations in the 4-year DMR maps will be published, together with some additional 
tests, in the June 10, 1996 issue of the Astrophysical Journal Letters. Copies of both of these 
papers are attached as an appendix to this report. 
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2. 3-point Correlations in the COBE-DMR Anisotropy Maps 

We have computed the 3-point correlation function of the COBE-DMR 2-year and 4-year 
anisotropy maps to search for evidence of skewness or asymmetry in the temperature fluctuations. 
Such behavior would have very significant implications for our understanding of the processes of 
galaxy formation, because our current models of galaxy formation predict that significant skewness 
should not be present in the DMR maps. 

Our results showed that the 3-point correlation function is consistent with zero and that 
the fluctuations are consistent with being Gaussian distributed. Moreover, with the improved 
sensitivity in the 2- and 4-year maps we were able to place significantly tighter upper limits 
on non-Gaussian fluctuations than previously existed, by nearly an order of magnitude. Our 
computation and analysis of 3-point correlations in the 2-year DMR maps was published in the 
Astrophysical Journal Letters, volume 446, page L67, 1995. Our computation and analysis of 
3-point correlations in the 4-year DMR maps will be published, together with some additional 
tests, in the June 10, 1996 issue of the Astrophysical Journal Letters. Copies of both of these 
papers are attached as an appendix to this report. 


3. 4-point Correlations in the COBE-DMR Anisotropy Maps 


We have computed a special case of the 4-point correlation function of the COBE-DMR 2-year 
maps to search for evidence of non-zero kurtosis in the temperature fluctuations: i.e. evidence that 
the “tails” of the temperature distribution differ significantly from a Gaussian distribution. Our 
results were inconclusive because of limitations with the DMR data: the level of “cosmic variance” 
is too high to make this higher-order statistic very useful. In other words, because of COBE’s 
coarse angular resolution, there are not enough independent points in the sky to usefully probe the 
tails of the temperature distribution. Since the results of this effort were not deemed interesting 
enough to submit for publication, we summarize the results of our study in some detail here. 


The simplest configuration of the 4-point function is one in which the 4 legs of the function 
are collapsed to two. This configuration depends only on the angular separation between the two 
legs, and is defined as 
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where t, is the temperature in pixel i, and the sum on j is restricted to pixels within angular 
separation bin a of pixel i. The ensemble average over realizations of the sky temperature is given 
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by 
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If the temperatures have a Gaussian distribution we may write 

(«?<?> = (<?)(<?> + 2(Wj) J 

= (C 2 (0)) 2 + 2(C 2 (o)) J 

where (C 2 (a)) is the ensemble-averaged 2-point correlation function. It follows that 

<C 4 (a)> = (C 2 (0)) 2 + 2 <C 2 (a)) 2 

The above result suggests that we define a generalized kurtosis statistic as 

A'(Q) = C 4 (a)-C 2 (0) 2 -2C 2 (o) 2 (2) 

One might expect this statistic to have zero mean for Gaussian distributed temperatures but 
that is not the case. The ensemble average of K is given by 

(K(a)) = (C 4 (a)) — (C 2 (0) 2 ) — 2(C 2 (a) 2 ) 

= (C 2 (0)) 2 + 2(C 2 (q)) 2 - <C 2 (0) 2 ) - 2 (C 2 (a) 2 ) 

which is, in general, not zero because of fluctuations in the 2-point function. If we define 

C 2 (a) = (C 2 (a)) + 6(a) 

where 5(a) is the deviation of an individual correlation function from the ensemble average, due, 
in general to cosmic variance and instrument noise, then, using (£(a)) ^ 0, we can write 

</v (a)> = (C 2 ( 0)> 2 + 2<c 2 (a)> 2 - <({C 2 (0)) + S( 0)) 2 > - 2{((C 2 (a)> + ^(a)) 2 ) 

= -(*(0) 2 )-2<*(a) 2 > 

Thus the larger the fluctuations in the 2-point function the more negative the mean of K will be. 
Because of this non- vanishing mean, it is problematic to interpret the results of a computation of 
K(a). In principle, one could avoid this by defining a modified kurtosis to be: 

K'(a) = C 4 («) - (C 2 (0)) 2 - (2 C 2 (a)) 2 (3) 


where the quantity subtracted off is a power of the ensemble- averaged 2-point function. This 
would have zero mean, but it requires that we know the ensemble-averaged 2-point function, a 
priori, which we don't. 
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Figure 1 shows the results of a Monte Carlo simulation of the 4-point function and the 
generalized kurtosis, as defined above. We simulated 2000 synthetic COBE-DMR sky maps with 
a Gaussian, scale-invariant, power-law model of CMB anisotropy, and instrument noise. For each 
realization, we computed the 4-point function, defined in (1), the generalized kurtosis, defined in 
(2), and the modified generalized kurtosis, defined in (3). Given 2000 realizations of each function 
we then compute its mean and standard deviation and plot the results as a function of angular 
separation. The solid lines in the plot show the mean and standard deviation of the 4-point 
function, C 4 (a). Note that the mean follows (C 2 (0) + 2C' 2 (a)), as expected, and has a rather large 
standard deviation, due mostly to cosmic variance. The short dashed lines show the mean and 
standard deviation of the kurtosis, K(a). Note that the mean is non-zero, as expected, but that 
the standard deviation is much smaller than with C 4 (a). Finally, the long dashed lines show the 
mean and standard deviation of the modified kurtosis, Ii '(a). Here the mean is zero (we have the 
advantage of knowing (C 2 (a)) because this is a simulation), but the standard deviation is much 
larger than with K(a). 

Even with the most sensitive statistic, K(a), the standard deviation is so large, ~ 30 /xK 4 , 
that it is completely insensitive to reasonable deviations from Gaussian statistics. To verify this, 
we have also run some simulations with “toy” non-Gaussian models in which the distribution of 
spherical harmonic coefficients, a( m , were chosen to be very different from Gaussian to see if the 
4-point statistic could pick this up. In general, the statistic was very insensitive to changes in the 
distribution of the a^ m ’s, further mitigating the interest of this statistic. 
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